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Abstract  --  A  new  data- treatment  procedure  allows  for  more  accurate 
determination  of  electron  temperatures  and  electron  concentrations  in 
analytical  plasmas.  A  Thomson- scattering  spectrum,  useful  for  these 
determinations,  is  often  not  purely  Gaussian  in  shape,  even  when  the 
probed  electrons  possess  a  Maxwellian  velocity  distribution. 

Nonetheless,  an  unambiguous  relationship  exists  between  electron 
temperatures  and  concentrations  that  truly  exist  in  a  plasma  and  those 
calculated  from  a  distorted  Thomson-scattering  spectrum.  Understanding 
this  relationship  permits  a  look-up  table  to  be  constructed,  from  which 
observed  values  can  be  corrected.  Theory  concerning  this  procedure  is 
described  and  details  tor  using  it  with  both  a  ruby  laser  and  frequency- 
doubled  Nd.YAG  laser  are  discussed.  Examples  of  electron  temperature 
and  electron  concentration  determined  with  this  procedure  in  an  ICP  are 
given.  The  possibility  of  using  the  new  procedure  to  study  electron- 
energy  distributions  is  also  assessed. 


% 
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1 .  INTRODUCTION 

The  inductively  coupled  plasma  (ICP)  has  become  an  extremely  important 
source  for  analytical  atomic  emission  and  mass  spectrometry  [1].  In  such  a 
plasma,  electrons  are  probably  the  most  active  and  important  species,  at 
least  as  far  as  energy  transfer,  excitation  and  ionization  are  concerned. 

It  is  principally  electrons  that  pick  up  energy  from  the  rf  field  generated 
by  the  induction  coil  and  transfer  it  to  ions  and  atoms  [2].  In  this  way, 
electrons  play  a  prominent  role  in  sustaining  the  plasma  and  in  transporting 
energy  throughout  the  discharge  [3,4].  As  importantly,  electrons  are 
involved  in  most  of  the  events  that  are  considered  significant  in  excitation 
and  ionization  of  analyte  species  [5-9] .  It  is  therefore  not  surprising 
that  electrons  have  been  studied  extensively  in  the  ICP  [10-24], 

Among  the  methods  to  investigate  electrons,  Thomson  scattering  [25] 
offers  a  number  of  advantages.  It  provides  values  for  both  electron 
concentration  and  temperature  simultaneously,  and  can  do  so  on  a  spatially 
and  temporally  resolved  basis.  Moreover,  under  proper  conditions,  it  can 
indicate  true  electron-energy  distributions.  Thomson  scattering  has  been 
used  as  a  diagnostic  tool  in  nuclear- fusion  research  for  more  than  20  years 

[26.27] .  More  recently,  the  feasibility  and  experimental  considerations  for 
using  Thomson  scattering  with  analytical  plasmas  were  reviewed  and  discussed 

[25.28] .  Preliminary  results  from  a  single-channel  scanning  device  and  a 
multichannel  instrument  have  proven  that  Thomson  scattering  can  be  a 
powerful  tool  for  probing  electron  temperatures  or  energy  distributions,  and 
electron  concentrations  in  the  ICP  [24,29,30]. 
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Unfortunately,  interpretation  of  Thomson- scattering  measurements  from 
an  ICP  is  more  complicated  than  those  in  a  fusion  plasma.  Fusion  plasmas 
normally  have  rather  high  electron  temperatures  (lO'^-lO^  eV)  and  relatively 
low  electron  number  densities  (10  -10  cm  ) ,  so  they  produce  a  Gaussian- 

shaped  Thomson- scattering  spectrum,  centered  at  the  incident  laser 
wavelength,  when  a  Maxwellian  velocity  distribution  of  electrons  exists 
[26,27,31].  Consequently,  a  plot  of  the  natural  logarithm  of  the  Thomson- 
scattering  signal  versus  the  square  of  the  wavelength  shift  produces  a 
straight  line,  from  the  slope  of  which  the  electron  temperature  can  be 
calculated  [26].  A  least-squares  linear  fit  is  usually  employed  in  this 
calculation,  since  experimental  errors  usually  occur.  In  addition,  the 
total  intensity  of  Thomson- scattered  light  can  be  used  to  determine  the 
value  for  a  localized  electron  number  density. 

Compared  to  fusion  plasmas,  analytical  plasmas  have  rather  low  electron 
temperatures  (<  1  eV)  and  relatively  high  electron  concentrations  (10^^ -10^^ 
cm  ).  Consequently,  a  Thomson-scattering  spectrum  from  an  ICP 
theoretically  will  not  have  a  purely  Gaussian  shape  even  if  the  electrons  in 
the  ICP  possess  a  Maxwellian  velocity  distribution.  A  large  observation 
angle  (approaching  backscattering)  and  a  short  incident  laser  wavelength 
(for  example,  a  frequency-doubled  Nd:YAG  laser  at  532  nm)  will  make  the 
scattering  spectrvim  approximate  a  Gaussian,  but  distortion  will  not  be 
completely  avoided.  One  can  still  apply  a  least-squares  fit  to  a  linearized 
version  of  such  a  spectrum  to  calculate  an  electron  temperature  and  can  use 
the  integrated  intensity  of  the  entire  spectrum  to  calculate  an  electron 
concentration  [24,25,29,30].  However,  the  distortion  of  the  spectrum  from  a 
true  C.i'issian  will  then  yield  an  overestimated  electron  temperature  and  an 
underestimated  electron  number  density;  the  scattering  signals  in  the  center 
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of  the  spectrum  are  depressed,  and  lead  to  an  erroneously  low  slope  and  a 
lower  overall  intensity  of  the  spectrum. 

In  the  present  paper,  an  improved  procedure  will  be  described  which 
serves  to  extract  valid  values  for  electron  temperature  and  electron  number 
density  from  an  experimental  Thomson- scattering  spectrum.  The  deviation 
between  these  values  and  those  determined  in  the  conventional  way  has  been 
found  to  depend  on  the  incident  laser  wavelength,  on  the  number  of  spectral 
intervals  (channels)  used  in  the  calculation,  and  on  the  magnitudes  of  the 
true  values  themselves.  Specifically,  a  longer  wavelength  of  the  incident 
laser,  a  smaller  number  of  spectral  channels,  a  higher  electron  number 
density  and  a  lower  electron  temperature  all  lead  to  a  more  substantial 
deviation  of  conventionally  determined  electron  temperatures  and 
concentrations  from  their  real  values.  These  points  will  be  illustrated  by 
means  of  data  obtained  from  recent  Thomson- scattering  measurements  on  an 
analytical  ICP.  Finally,  the  feasibility  of  using  the  new  procedure  for 
detecting  a  non -Maxwellian  electron  velocity  distribution  will  be  assessed. 

2.  THEORETICAL  CONSIDERATIONS 

When  an  intense  light  beam  strikes  a  plasma,  free  electrons  are 
accelerated  in  resonance  with  the  incident  electromagnetic  field.  As  a 
consequence,  the  oscillating  electrons  act  as  miniature  sources  themselves 
and  emit  radiation  of  their  own.  This  secondary  radiation  is  termed  Thomson 
scattering.  When  the  electrons  are  in  rapid  motion,  as  is  certainly  the 
case  in  a  hot,  rare-gas  plasma,  the  scattered  radiation  will  be  Doppler- 
shifted  and  will  produce  a  so-called  Thomson- scattering  spectrum,  from  which 
electron  energies  and  concentrations  can  be  deduced  [24,25,29], 


6 


Of  course,  ions  experience  the  same  electromagnetic  field  as  do 
electrons.  However,  their  mass  is  so  large  that  they  are  displaced  very 
little  and  therefore  scatter  relatively  weakly.  Nonetheless,  ions  can  under 
some  conditions  greatly  influence  the  electron  density  oscillation,  and 
thereby  lead  to  a  change  not  only  in  the  shape  of  the  Thomson- scattering 
spectrum  but  also  in  the  total  intensity  of  scattered  light.  This  influence 
makes  interpretation  of  a  Thomson  spectrum  difficult  and  often  leads  to 
errors  in  calculated  electron  temperature  and  number-density  values. 

To  see  how  these  errors  can  be  minimized,  let  us  briefly  review  the 
quantitative  aspects  of  Thomson  scattering.  The  spatial  and  frequency 
(wavelength)  distribution  of  Thomson- scattered  light  is  expressed  by  the 
differential  scattering  cross  section  a  (k,  w)  [26,32]: 

C7(k,  w)  =  <7.j.‘S(k,  w) ,  (1) 

where 

=  (e^/m^c^)^  sin^  7  (2) 

is  the  Thomson  cross  section,  S(k,  w)  the  so-called  dynamic  form  factor,  k 
the  differential  scattering  wave  vector  [2.5,26],  and  u>  the  frequency  shift 
of  the  scattered  light  from  the  incident- light  frequency.  In  Eqn.  (2),  e, 
and  c  are  the  charge  and  mass  of  an  electron  and  the  speed  of  light, 
respectively,  and  7  is  the  angle  between  the  observation  direction  and  the 
electric  vector  of  the  incident  light  beam.  The  maximum  scattering  signal 
is  produced  at  an  angle  equal  to  90°. 


For  a  Maxwellian  electron- velocity  distribution,  the  dynamic  form 
factor  can  be  decomposed  into  two  terms,  the  first  of  which  deals  with  free 


electron  scattering  and  the  second  of  which  describes  the  effect  of 
neighboring  ions  on  electron  scattering  [26,32]: 


a2  2 

s(k,w)  dw  -  r^(x  )dx  +  z  ( - ~)^  r3(x.)dXj^  : 

l+a^ 


(3) 


ra(Xe) 


exp(-Xg‘^) 


w(x  )]• 


;  (Xi) 


expC-X'^j^) 


y;^|i+^2„(Xj^) 


(^) 


and 


0) 

Xe  =  -  ;  X.  =  - 


kV, 


kVi 


(5) 


,v2 

o  a  e 

P-'Z- - ^  •  - 

1+Q^  T,. 


(6) 


W(X)  =  1  -2Xe 


■X 


x  2 
P 


e*^  dp  -  ijlr  Xe' 


(7) 


kAr 


(8) 


is  the  scattering  parameter.  Note  from  the  shape  of  (X^) 

that  S(k,  w)  is  the  sum  of  two  Gaussians  when  a  «  1.  In  Eqns  (3)  and 
(6),  Z  is  the  charge  on  the  ion.  In  Eqn  (5)  X^  and  Xj^  are  dimensionless 
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frequencies,  and  and  are  the  most  probable  velocities  of  the 
electron  and  ion,  respectively,  which  can  be  determined  from 


where  ko  is  the  Boltzmann  constant,  T,,  and  T,-  are  the  electron  and  ion 
rS  ’  e  1 

temperatures,  respectively,  and  m^  is  the  mass  of  the  ion.  In  Eqn  (8), 
Ajj  is  the  Debye  length  given  by 


(10) 


where  n^  is  the  electron  number  density.  The  Debye  length  Aj^  is  a 
measure  of  the  charge -shielding  distance  in  a  plasma  [33].  Electrons 
exhibit  individual  behavior  on  a  scale  much  shorter  than  the  Debye 
length  and  can  be  considered  free;  collective  electron  motion  occurs 
over  a  range  longer  than  the  Debye  length,  and  results  in  correlated  or 
coherent  phenomena. 

The  scalar  quantity  k  in  Eqns  (5)  and  (8)  has  to  do  with  the 


distance  in  the  plasma  over  which  scattering  occurs.  When  a  laser  beam 
of  wavelength  A^  is  incident  on  a  plasma,  light  scattering  at  an  angle  8 
with  respect  to  the  incident  direction  occurs  on  a  distance  scale  of 
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1/k,  where  k  is  the  magnitude  of  the  vector  quantity  k  and  is  given  by 
[26,32] 

4jrsin(5/2) 

k - -  .  (11) 

■'o 

Clearly,  when  (l/k)«Xjj,  a«l  (see  Eqn  (8)),  and  electrons  in  the 

plasma  would  behave  as  if  they  were  free  in  the  electromagnetic  field  of 

the  laser  beam.  As  a  result,  the  spectrum  of  scattered  light  would 

directly  reflect  the  electron-velocity  distribution.  This  behavior  can 

be  readily  seen  from  Eqns  (3)  and  (4),  because  the  second  term  on  the 

right  hand  side  of  Eqn  (3)  can  then  be  neglected,  as  can  the  value  of 
2 

a  W(Xg)  in  the  first  term.  Under  these  conditions, 

exp(-Xg^) 

S(k,w)  <k)  -  -  dX  (12) 

R 

which  shows  that  the  scattering  spectrum  will  have  a  Gaussian  shape. 

Integrating  Eqn  (12)  over  frequency  w  for  the  entire  spectrum  leads 
to 


Combining  Eqns  (1)  and  (13)  indicates  that  the  integral  of  the 
differential  scattering  cross  section  a  (k,cj)  over  frequency  becomes 
equal  to  the  Thomson  cross  section,  when  a«l .  Therefore,  the  detected 
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scattering  signal  I(w)  is  proportional  to  S(k,a>)  given  by  Eqn  (12).  From 
this  relationship  and  Eqn  (12) ,  one  can  obtain  the  following  expression 
by  taking  a  logarithm: 


2 

From  Eqn  (15),  the  slope  of  a  plot  of  In  1(A)  vs.  (AA)  can  be  used 
to  determine  the  electron  temperature  T^.  This  is  the  usual  way  to  find 
the  electron  temperature  from  a  Thomson- scatter ing  spectrum,  and  is 
convenient  and  accurate  when  a  «  1.  Furthermore,  under  these 
conditions  (q«1)  the  calibrated  signal  can  be  used  to  determine  the 
electron  concentration  precisely,  because  the  average  light-scattering 
cross  section  of  an  electron  is  simply  equal  to  the  Thomson  cross 
section  as  discussed  above.  As  a  result,  the  absolute  scattering 
intensity  is  directly  proportional  to  the  electron  number  density. 
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Errors  Resulting  from  a  /  0 

Unfortunately,  the  condition  a«l  does  not  exist  in  the  ICP.  For  a 
given  electron  temperature  and  electron  concentration  n^,  a  shorter 
laser  wavelength  and  larger  observation  angle  6  will  produce  a 
smaller  a  value  (cf.  Eqns  (8)  and  (11)).  However,  even  for  a 
backscattering  geometry  (6  =  180°),  the  wavelength  of  a  ruby  laser 
(694  nm)  or  a  frequency  doubled  Nd:YAG  laser  (532  nm)  will  yield  a  values 
that  lie  between  0.1  and  1.  Consequently,  the  electrons  will  not  behave 
independently  and  the  scattering  spectra  will  deviate  from  their  ideally 
Gaussian  shape. 

Electron-Temperature  Errors 

The  degree  of  this  deviation  and  its  dependence  on  a  are  revealed 
in  the  curves  of  Fig.  1  calculated  according  to  Eqns  (4)  and  (7). 
Interestingly,  the  curves  all  nearly  cross  at  about  =  1.0.  At  values 
of  Xg  below  1.0,  the  center  of  each  curve  is  depressed  when  a>0  and  by 
an  amount  that  increases  with  a.  In  contrast,  the  tail  (X^  >  1.0)  of 
each  spectrum  is  raised  slightly  when  a  is  above  zero. 

Because  only  the  curve  in  Fig.  1  corresponding  to  o  =  0  is 

2 

Gaussian,  plotting  of  In  r(;j(Xg)  vs.  Xg  for  the  other  curves  results  in 
bending  (see  Figs.  2  and  3).  To  construct  Figs.  2  and  3,  ten  values  of 
Xg  ranging  from  0.2  to  2.0  (separated  by  0.2)  were  used,  and  correspond 
roughly  to  the  channels  used  in  a  multichannel  Thomson- scatter ing 
experiment.  At  each  of  these  values,  the  quantity  In  rjj(Xg)  was 
calculated  and  plotted  as  circles  (Fig.  2)  or  triangles  (Fig.  3) .  A 
least  -  squares  linear  fit  to  these  plotted  points  (solid  line)  can  then 
tie  compared  to  the  straight  line  that  would  have  been  produced  by  a  =  0 
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(dashed  line) .  It  can  be  seen  that  the  slope  of  the  least  -  squares 
linear  fit  in  both  Figs.  2  and  3  (solid  line)  is  less  steep  than  that  of 
the  plot  for  a  =  0  (dashed  line).  As  a  consequence,  the  electron 
temperature  calculated  from  Eqn  (14)  would  be  overestimated.  Moreover, 
the  larger  the  a  value,  the  greater  would  be  the  error  in  computed 
electron  temperature  (compare  Figs.  2  and  3). 

Note  that  the  disparity  among  the  curves  in  Fig.  1  is  relatively 
smaller  in  the  tail  than  in  the  central  part.  This  relationship  is 
clear  also  in  Figs.  2  and  3,  especially  if  one  examines  the  deviations 
of  circles  or  triangles  from  the  dashed  lines.  Obviously,  using  only 
scattering  signals  from  the  tail  channels  of  a  Thomson- scattering 
spectrum  results  in  a  least-squares-linear-fit  value  for  electron 
temperature  that  is  closer  to  the  real  electron  temperature  than  that 
obtained  when  the  entire  scattering  spectrum  is  employed  [28]. 
Unfortunately,  the  scattering  intensity  decreases  so  quickly  with 
wavelength  shift  (see  Fig.  1)  that  the  detected  signals  in  far-wing 
channels  are  usually  very  noisy.  It  is  clearly  undesirable  not  to  use 
scattering  signals  from  channels  located  closer  to  the  incident  laser 
wavelength . 

Errors  in  Electron  Number  Density 

When  the  scattering  parameter  a  is  not  nearly  zero,  as  is  the  case 
in  the  ICP,  the  average  light-scattering  cross  section  of  an  electron  is 
not  equal  to  the  Thomson  cross  section,  but  rather  decreases  with  the  a 
value  itself.  From  Eqns  (1),  (3),  and  (4),  one  notes  that  the 
scattering  cross  section  of  an  electron  when  q  >  0  is  smaller  than  the 
Thomson  cross  section  by  a  factor  of  [25,26] 
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a(k,w)dw 


*qd 

S(k,6i))  da>  ■= 

J  _00 


(Xe)  dX^ 


dX^  1 

„  y^|l+a2y(Xg)|2  l+a2 


(16) 


This  fraction  is  therefore  equal  to  the  area  between  the  axis  and  the 
curve  with  a  specified  a  value  in  Fig.  1,  over  the  entire  range  of  . 
Clearly,  the  larger  the  a  value,  the  smaller  the  area. 

Because  of  this  fixed  relationship,  one  could  calculate  valid 

2 

numbers  for  the  electron  number  density,  just  by  using  o  —  o^/(l  +  a  ) 
instead  of  for  the  average  scattering  cross  section  of  the  electron. 
The  difficulty  encountered,  however,  is  that  the  correct  a  value  is 
unknown  before  a  real  electron  temperature  and  electron  concentration 
have  been  determined;  in  turn,  the  determination  of  an  electron 
temperature  and  concentration  requires  knowledge  of  the  a  value. 


Correction  Procedure 

To  overcome  this  cyclic  and  apparently  insoluble  problem,  we  must 
recognize  that  there  is  an  unambiguous  relationship  among  a  specific 
electron  temperature,  an  electron  concentration,  and  the  a  value  [see 
Eqns  (8),  (10),  and  (11)].  Moreover,  all  affect  a  Thomson- scattering 
spectrum  in  well  defined  ways  (see  Eqns.  (4),  (5),  and  (7)]. 
Consequently,  for  every  combination  of  T^  and  n^ ,  there  will  be  a 
particular  a  and  a  specific  Thomson- spectrum  shape. 
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Practically,  one  can  proceed  by  assuming  values  for  and  n^  and 
by  calculating  from  them  a  number  for  a.  These  parameters  then  permit  a 
complete  Thomson- scattering  spectrum  to  be  computed  and,  from  it,  a 
least-squares  linear  fit  and  an  integrated  spectrum.  The  resulting 
slope  and  integral  can  be  used  as  formerly  to  derive  "experimental" 
representations  for  T^  and  n^.  This  computational  procedure  can  be 
repeated  for  all  anticipated  or  realistic  combinations  of  electron 
temperatures  and  number  densities  expected  in  the  ICP,  to  yield  a  matrix 
or  table  that  can  be  employed  to  equate  "experimental"  and  real  values. 
The  table  prepared  in  this  way  can  then  be  used  in  reverse,  to  deduce 
true  or  correct  values  from  those  originally  determined. 

3.  PREPARATION  AND  USE  OF  CORRECTION  TABLES 
Correction  tables  such  as  described  above  can  be  computed  for 
specific  laser  wavelengths,  adopted  observation  angles,  6,  and 

particular  numbers  of  detection  channels,  each  of  which  covers  a 
specified  wavelength  range.  In  this  way,  a  correction  table  can  be 
tailored  to  an  individual  experimental  arrangement.  Table  1  and  Table  2 
are  two  examples  that  have  proven  useful  in  our  own  studies.  Table  1 
has  been  computed  for  a  ruby  laser  excitation  source  and  with  each 
detection  channel  covering  0.45  nm;  Table  2  applies  to  the  use  of  a 
frequency-doubled  Nd:YAG  laser,  with  each  spectral  channel  covering 
0.34  nm.  Both  tables  pertain  to  an  observation  angle  of  135”  and  assume 
that  nine  spectral  channels  are  used.  Although  these  tables  have  proven 
useful  in  practice,  they  can,  of  course,  be  prepared  with  finer 
resolution  in  both  electron  temperature  and  number  density.  In 
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situations  where  signal-to-noise  ratio  is  extremely  high,  such  smaller 
increments  would  obviate  the  need  for  interpolation. 

Effect  of  Number  of  Spectral  Channels 

As  indicated  earlier,  the  electron  temperature  deduced  from  a 
Thomson- scattering  spectrum  is  influenced  strongly  by  intensity  values 
on  the  "wings"  of  the  spectrum,  whereas  a  calculated  electron 
concentration  is  affected  mostly  by  values  near  the  central  maximum  of 
the  spectruun.  It  is  therefore  not  surprising  that  the  number  of 
spectral  channels  used  in  a  determination  of  and  n^  would  affect  the 
two  differently,  assuming  each  channel  covers  the  same  spectral 
interval.  This  behavior  is  explored  in  Figs.  4-7.  In  these  Figures, 
the  number  of  channels  that  is  cited  refers  to  those  in  the  center  of 
the  spectrum;  however,  they  begin  just  beyond  the  central  (Rayleigh) 
channel,  which  has  been  left  vacant. 

From  Figs.  4  and  5,  it  is  apparent  that  the  deviation  of  an 
observed  T^  from  the  real  T^  decreases  rapidly  as  more  channels  are 
used.  Also,  because  it  produces  a  narrower  Thorason-scattering  spectrum, 
the  frequency-doubled  Nd;YAG  laser  produces  a  smaller  deviation  than  the 
ruby  laser  for  any  particular  number  of  channels.  Not  surprisingly,  the 
deviation  is  very  sensitive  to  electron  concentration  and  is  greatest  at 
large  Og. 

Figs.  6  and  7  show  "experimental"  values  of  n^  calculated  as  a 
function  of  the  number  of  channels  used  for  the  same  sets  of  real  T^  and 
n^  as  in  Figs.  4  and  5.  Not  surprisingly,  the  disparity  between 
observed  and  real  n^  values  remains  almost  unchanged  when  more  than  7 
channels  are  used;  the  signal  from  the  tail  channels  of  the  spectrum  i.s 


so  small  that  it  does  not  contribute  appreciably  to  the  calculated 
electron  concentration  (see  Fig.  1).  Moreover,  both  the  ruby  and 
frequency -doubled  Nd.YAG  lasers  provide  about  the  same  accuracy  for  n^ 
determination,  particularly  when  n^  is  low  (see  Fig.  7).  Importantly, 
however,  the  error  in  deduced  n^  (Figs.  6  and  7)  is  never  as  low  as  that 
in  Tg  (Figs.  4  and  5),  regardless  of  how  many  observation  channels  are 
employed.  This  point  will  be  discussed  later. 

Application  of  Correction  Tables 

The  above-mentioned  procedure  for  tracing  real  electron 
temperatures  and  concentrations  has  been  applied  to  some  of  our  recent 
Thomson- scatter ing  data  and  has  proven  effective  and  useful.  Examples 
are  given  in  Figs.  8  and  9.  Experimental  conditions  used  to  obtain 
these  data  are  cited  in  Table  3.  No  sample  aerosol  was  used.  In  Figs. 

8  and  9,  the  open  circles  represent  values  calculated  simply  by  using  a 
least- squares  linear  fit  for  T^  and  the  Thomson  cross  section  for  n^ . 

The  filled  circles  are  the  corresponding  values  traced  using  Table  1  (or 
similar)  and  should  indicate  the  true  electron  temperatures  and  number 
densities.  In  both  figures,  nine  spectral  channels  were  used  for  data 
obtained  at  observation  heights  of  5.0.  7.5,  and  10.0  mm  above  the  load 
coil  (ALC) ;  seven  and  six  channels  were  used  for  data  taken  at  12.5  and 
15.0  mm  ALC,  respectively.  Fewer  channels  were  employed  higher  in  the 

plasma  because  of  the  lower  n^  that  exists  there;  this  lower  n^  yields  a 

weaker  Thomson  signal  and  produces  a  poorer  s ignal - to-noise  ratio  in  the 
tail  channels  of  the  spectrum. 

Figures  8  and  9  indicate  that  the  new  correction  procedure  becomes 

more  important  for  both  T^  and  n^  at  lower  observation  heights.  The 
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reason  is  that  the  electron  number  density  is  quite  high  at  positions 
low  in  the  plasma  and  produces  a  larger  a  value.  In  turn,  a  large  a 
leads  to  a  more  serious  depression  of  the  central  channel  signals  and  a 
greater  distortion  of  the  Thomson  spectrum.  Interestingly,  the 
corrected  data  in  Figs.  8  and  9  indicate  that  the  electron  temperature 
does  not  drop  as  rapidly  with  observation  height  as  does  electron  number 
density.  It  would  seem  then  that  electron  temperature  is  not 
responsible  for  causing  major  changes  in  the  a  value  in  different 
regions  of  the  ICP.  In  general,  the  corrected  data  in  Figs.  8  and  9  are 
similar  to  those  previously  reported  (12,13,16,18,19,22).  Although  the 
electron  temperatures  cited  in  Fig.  8  are  about  20%  lower  than  those  by 
Batal  et  al.  [16],  their  observation  height  was  below  ours  (2  mm  above 
the  load  coil)  and  their  input  RF  power  was  1.5  kW  (compared  to  our 
1  kW) .  Also,  the  electron  number  densities  in  Fig.  9  are  roughly  twice 
those  reported  by  Caughlin  and  Blades  (18)  at  low  positions  in  the 
plasma  (5  mm  above  the  load  coil),  but  are  in  good  agreement  at  higher 
spatial  points  (15  mm  above  the  load  coil).  This  disparity  might  be  due 
to  the  slightly  different  torches  used  in  the  two  studies. 

It  should  be  pointed  out  that  Caughlin  and  Blades  used  the  Vidal, 

Cooper  and  Smith  tabulations  of  half-width  to  deduct  their  n^  values;  if 
they  had  used  Goode  and  Deavor's  Stark  algrithm  [34],  their  n^  would  be 
lower  than  ours  even  at  an  observation  height  of  15  mm  above  the  load  coil. 

From  the  corrected  data  shown  in  Figs.  8  and  9,  it  is  clear  that  the 
electron  temperatures  are  not  as  high  as  were  previously  thought.  In  other 
words,  the  ICP  is  probably  not  very  far  from  LTE,  and  therefore  the  close- 
to-LTE  concept  is  probably  useful  in  evaluating  ICP  behavior  [5,35,36).  In 
fact,  our  Tp  values  are  slightly  lower  than  the  LTE  values  determined  by  our 
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measured  ,  obtained  using  the  Saha  equation.  This  relationship  reflects 
the  recombining  character  of  the  plasma,  and  would  lead  to  lower  ion/atom 
emission- intensity  ratios  than  only  n^-determined  LTE  conditions  would 
predict,  as  has  been  observed  by  Caughlin  and  Blades  [36], 

4.  DISCUSSION  AND  CONCLUSIONS 

When  electron  number  density  is  high,  the  deviation  of 
conventionally  calculated  data  from  real  values  is  large  for  both  and 
n^  (see  Figs.  4,  6,  8,  and  9);  the  depression  of  the  central  channels  in 
a  Thomson  spectrum  induced  by  a  large  a  value  clearly  influences  both  T^ 
and  n^  calculations  greatly.  In  contrast,  when  the  electron  number 
density  is  low,  the  deviation  in  calculated  T^  is  greatly  reduced,  as 
can  be  seen  by  comparing  Figs.  4  and  5  (see  also  Fig.  8).  However,  the 
deviation  in  n^  is  high  even  at  low  n^  (compare  Figs.  6  and  7,  and  see 
also  Fig.  9).  The  reason  is  that  the  underestimation  of  n^  is  due  not 
only  to  the  a  value,  but  also  to  the  fact  that  the  contribution  of  the 
central  (Rayleigh)  channel  to  the  Thomson  spectrum  cannot  be  accurately 
assessed.  At  that  wavelength,  Rayleigh  scattering  and  stray  light  are 
normally  so  high  that  their  contribution  to  the  observed  signal  cannot 
be  subtracted.  When  electron  nvimber  density  is  low,  the  effect  of  a 
Large  a  value  on  an  n^  calculation  becomes  negligible,  but  the  central - 
channel  signal  becomes  more  important.  In  turn,  ignoring  the  central - 
channel  signal  leads  again  to  an  underestimation  of  n^ .  Importantly, 
this  error  is  automatically  avoided  when  the  new  correction  procedure  is 
performed . 

Although  the  experimental  precision  of  our  current  Thomson  scattering 
measurements  is,  on  the  a'-'erage,  around  10%,  we  feel  that  the  procedure  for 
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tracing  a  true  and  n^  described  above  should  be  applicable  to  Thomson- 
scattering  data  obtained  in  most  regions  of  analytical  interest  in  the  ICP. 
Nonetheless,  the  accuracy  of  the  resulting  values  still  depends  on  how  well 
the  relative  sensitivities  of  different  spectral  channels  have  been 
normalized  and  how  well  the  absolute  intensity  of  the  entire  spectrum  has 
been  calibrated.  For  data  obtained  on  a  well  normalized  and  calibrated 
Thomson-scattering  apparatus  with  high  signal -to-noise  ratio,  the  new 
procedure  should  offer  faithful  values  of  the  true  and  n^ .  Of  course , 
the  new  procedure  might  not  be  the  only  way  to  trace  a  real  Tg  and  n^ ;  curve 
fitting  by  a  computer  is  an  alternative  and  attractive  approach  that  we  are 
investigating. 

Another  approach  exists  for  the  unambiguous  measurement  of  real  Tg  and 
Og  by  Thomson  scattering,  taking  advantage  of  a  »  1  conditions  [26,28].  In 
this  method,  Og  is  measured  directly  by  obseirving  the  resonance  satellite 
related  to  the  electron  plasma  frequency,  using  the  relation 


0.2  =  (1  +  3/„2) 


(17) 


where  w  is  the  satellite  peak  frequency  and  Wp  is  the  electron  plasma 
frequency,  which  depends  only  upon  the  ng  value 


(18) 


At  very  shallow  scattering  angles,  a  »  1  and  w  =  Wp,  so  ng  can  be 

o 

calculated  directly  from  Eqn.  (18).  Then  at  a  steeper  angle  (say  30  ) ,  w 
can  be  measured  again,  permitting  a  to  be  computed  from  Eqn.  (17)  and  Tg  to 
be  determined  from  Eqns .  (8),  (10),  and  (11).  However,  this  approach 
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requires  a  very  small  observation  angle  6,  and  at  least  one  additional 
observation  direction  for  the  determination  of  both  and  n^. 

Consequently,  more  effective  stray  light  rejection  and  a  more  sophisticated 
experimental  arrangement  would  be  necessary  than  in  the  present  Thomson- 
scattering  system. 

Application  to  Non -Maxwellian  Electron  Energy  Distributions 

In  most  regions  of  the  ICP,  a  Thomson- scattering  spectrum  will  not 
have  a  purely  Gaussian  shape,  and  it  is  appropriate  to  ask  whether  such 
a  spectrum  can  be  used  to  study  electron-energy  distributions.  Indeed 
it  can,  but  only  if  the  scattering  system  has  been  well  normalized,  and 
offers  a  high  s ignal - to -noise  ratio. 

In  an  ICP,  where  the  value  of  q  departs  substantially  from  zero,  a 
departure  from  Maxwellian  behavior  could  be  discerned  in  the  following 
way.  First,  the  "experimental"  values  of  and  n^  can  be  calculated  in 
the  conventional  manner  and  the  corresponding  true  values  found  by 
performing  the  procedure  suggested  in  this  paper.  The  theoretical 
signal  magnitude  from  each  channel  used  in  the  Thomson- scatter ing 
measurements  can  then  be  calculated  by  using  the  true  (corrected)  T^  and 
n^  magnitudes  and  Eqns  (4)  and  (7).  If  the  electron  velocity 
distribution  is  Maxwellian,  the  calculated  signal  in  each  spectral 
channel  should  be  in  good  agreement  with  the  corresponding  measured  one. 
Otherwise,  some  of  the  calculated  signal  intensities  would  be  higher  and 
some  would  be  lower  than  the  corresponding  measured  data,  but  the  least- 
squares  linear  fit  should  be  almost  the  same  for  both  calculated  and 
measured  intensities.  Carefully  examining  how  the  measured  signals 
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deviate  from  the  calculated  intensities  would  give  information  about  the 
electron  energy  distribution  in  the  plasma. 

Finally,  it  should  be  pointed  out  that  the  procedure  suggested  here 
for  tracing  real  Tg  and  real  n^  from  observed  values  applies  not  only  to 
the  ICP,  but  also  to  any  other  plasmas  with  comparable  ranges  of  Tg  and 

’"e- 
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Table  1.  Correction  table  for  ruby-laser  source.  Real  (true)  (top  line) 

and  Ng  (left  column)  values  and  corresponding  observed  (experimental)  values 
(table  entries)  pertaining  to  use  of  a  ruby  laser  and  nine  0.45-nm-wide  spectral 
channels.  The  range  of  channels  cover  a  wavelength  shift  (at  the  center  of 
each  channel)  from  0.45  to  4.05  nm.  Deleted  channels  listed  on  the  upper-left 
corner  of  the  table  are  available  in  the  Thomson-scattering  system,  but  are  not 
used  for  reasons  described  in  the  text. 


channel  s  (ruby)  :  0  10  11  12  13 


fftAL  Uo 

REAL 

Te 

0.6&OO0O£*04 

0.60OOOO£-*04 

0.6500006*04 

e.reeeeaE.e'* 

0.7500006*04 

0.8000006*04 

e.eeeeeoE^e* 

0.9000006*04 

e.l00000£«16 

T«=  0.673268£*e4 

Ne=  0.767741£*16 

0,6273666*04 

0.7788006*16 

0.6816076*04 

0.7663476*15 

0.7359386*04 

0,7965986*15 

0.7903256*04 

0.8037226*15 

0.8447356*04 

0.6098566*15 

0.8991426*04 

0.8161136*15 

0.9535236*04 

0.8196896*15 

e.i6000C£«ie 

Te-  0.6ed39O£*04 
NE=  0.109S06e*16 

0.6427916*04 

0.1114636*16 
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Table  2.  Correction  table  for  Nd/YAG  laser.  Real  (true)  T^  (top  line)  and 
(left  column)  values  and  corresponding  observed  (experimental)  values  (table 
entries)  pertaining  to  use  of  a  frequency-doubled  NdrYAG  laser,  and  nine  0.34- 
nm-wide  spectral  channels,  covering  a  wavelength  shift  (at  the  center  of  each 
channel)  from  0.39  to  3.08  nm.  Deleted  channels  listed  on  the  upper  left  corner 
of  the  table  are  available  in  the  Thomson-scattering  system,  but  are  not  used 
for  reasons  described  in  the  text. 
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Table  3.  ICP  operating  conditions 


Torch 

RF  frequency 
RF  input  power 
Outer  argon  flow 
Intermediate  argon  flow 
Central  flow 
Sample  flow  rate 


Commercial  MAK  torch 
27.12  MHz 
1.0  kw 
10  1/min 
0.5  1/min 
0 . 5  1/min 
none 
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Figure  Captions 

Figure  1.  Calculated  function  r£,(Xg)  for  several  values  of  the 

scattering  parameter  a.  The  magnitude  of  is 

roughly  proportional  to  the  expected  spectral  intensity 
of  Thomson  scattering.  For  a  given  experimental  system 
and  scattering  geometry,  is  proportional  to  frequency 
shift  from  that  of  the  incident  laser  (see  Eqns  (3)  (4) , 
and  ( 7 ) ) . 

Figure  2.  Effect  of  non-Gaussian  Thomson- scattering  spectrum  on 

errors  produced  by  a  linear- fitting  process.  Calculated 
values  of  In  r^(Xg)  for  a  -  0.3;  filled  circles;  least- 
squares  linear  fit  of  the  calculated  points:  solid  line; 
calculated  function  In  ^^(Xg)  for  a  —  0:  dashed  line. 

Figure  3.  Effect  of  non-Gaussian  Thomson- scattering  spectrum  on 

errors  produced  by  a  linear- fitting  process.  Calculated 
values  of  In  r^(Xg)  for  a  -  0.5;  filled  triangles;  least- 
squares  linear  fit  of  the  calculated  points:  solid  line; 
calculated  function  In  ^'^((Xg)  for  q  —  0;  dashed  line. 

Figure  4.  "Experimental"  electron  temperatures  calculated  as  a 

function  of  the  number  of  spectral  channels  used:  open 
circles  for  ruby  laser;  filled  circles  for  frequency- 
doubled  Nd.YAG  laser.  A  real  Tg  -  8000  K  (dashed  line) 

1C 

and  a  real  N  =  3  x  10  cm  are  assumed.  Number  of 
e 
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Figure  5 . 


Figure  6. 


Figure  7. 


channels  corresponds  to  the  width  of  observed  spectrum: 
each  channel  is  0.45  nm  in  the  case  of  the  ruby  laser  and 
0.34  nm  for  the  Nd:YAG  laser. 

"Experimental"  electron  temperatures  calculated  as  a 
function  of  the  number  of  spectral  channels  used:  open 
circles  for  ruby  laser;  filled  circles  for  frequency- 
doubled  Nd:YAG  laser.  A  real  —  6000  K  (dashed  line) 
and  a  real  n^  —  lO'"'  cm*  are  assumed.  Number  of 
channels  corresponds  to  the  width  of  observed  spectrum: 
each  channel  is  0.45  nm  in  the  case  of  the  ruby  laser, 
and  0.34  nm  for  the  YAG  laser. 

"Experimental"  electron  concentrations  calculated  as  a 
function  of  the  number  of  spectral  channels  used;  open 
circles  for  ruby  laser;  filled  circles  for  frequency- 
doubled  Nd:YAG  laser.  A  real  “  8000  K  and  a  real  n^  = 
3  X  10  cm  (dashed  line)  are  assumed.  Number  of 
channels  corresponds  to  the  width  of  observed  spectrum; 
each  channel  is  0.45  nm  in  the  case  of  the  ruby  laser, 
and  0.34  nm  for  the  YAG  laser. 

"Experimental"  electron  concentrations  calculated  as  a 
function  of  the  number  of  spectral  channels  used:  open 
circles  for  both  ruby  laser  and  frequency-doubled  Nd:YAG 
laser.  A  real  =  6000  K  and  -  10^^  cm’^  (dashed 


line)  are  assumed.  Number  of  channels  corresponds  to  the 


width  of  observed  spectrum;  each  channel  is  0.45  nm  in 
the  case  of  the  ruby  laser,  and  0.34  nm  for  the  YAG 
laser. 


Figure  8 . 


Figure  9 . 


Electron  temperatures  measured  by  Thomson  scattering 
from  an  ICP  at  different  observation  heights  and  at  3  mm 
off  the  plasma  axis.  Open  circles:  observed  values; 
filled  circles;  traced  values  (see  text  for  detail). 

ICP  operating  power  -  1  kW;  no  aerosol  is  being 
introduced.  See  Table  3  for  other  operating  conditions. 

Electron  concentrations  n^  measured  by  Thomson  scattering 
from  an  ICP  at  different  observation  heights  and  at  3  mm 
off  the  plasma  axis.  Open  circles:  observed  values; 
filled  circles:  traced  values  (see  text  for  detail). 

ICP  operating  power  -  1  kW;  no  aerosol  is  being 
introduced.  See  Table  3  for  other  operating  conditions. 


Number  of  channels  used  in  calculation 


Number  of  channels  used  for  calculation 


1.2E15 


Number  of  channels  used  for  calculation 


Observation 


6.0E15 


5  5.0  7.5  10.0  12.5  15.0  i: 

Observation  height  above  the  load  coil  (mm) 


